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ABSTRACT 


Several small-scale, special-purpose computer programs that were 
written to provide design data for a pulse-slimming filter over a wide 
range of parameter variations are contained herein. A large-scale 
computer program is also used to provide additional data. Some syn- 
thesis techniques are discussed, and some of the computer program 
outputs are directly applied in illustrating synthesis techniques for 


a specific pulse-slimming transfer function. 
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I. INTRODUCTION 


Network analysis and synthesis, in general, continues to become less 
time-consuming and complex as a result of the continued upgrading of 
large-scale network-analysis computer programs such as those listed in 
Ref. 1 and the increased interest shown in computer-aided design by 
engineers at all levels. 

A major disadvantage of the large-scale network analysis programs, 
however, is the vast amount of core storage usually required by them. 

If a large computer is not available to the design engineer, many of the 
more sophisticated programs will not be at his disposal. Additionally, 
owing to the fact that network synthesis, particularly filter design, is 
approached from many aspects depending on the specific application, the 
available large-scale programs must frequently be supplemented by a host 
of smaller, special-purpose programs in order to provide the required 
design information. Thus, it is essential that the present-day design 
engineer be well versed in network analysis programs and computer capabil- 
ities in order to successfully compete with his contemporaries. 

Computer-aided design procedures are exemplified herein by a compre- 
hensive analysis of the pulse-slimming filter designed by H. M. Sierra 
[Ref. 2]. The use of small, special-purpose programs coupled with a 
large-scale program provides the required design parameters. 

some aspects of lattice synthesis are discussed, and with the aid of 
the computer programs some of these synthesis techniques are applied to 


the Sierra transfer functions. 


A. THE SIERRA FILTER 

Reference 2 contains a detailed description of the Sierra Filter; 
however, for completeness, a brief description of this filter will be 
included. 

Consider an input pulse, f.(t), approximated by the normal curve 


Z ey 


fe) = ere cent sa<t<a (ia 


0 otherwise. 


The error in this approximation at t = a is simply the value of the 
Gaussian function at t =a or 


€ = Sync : hes) 


The problem as outlined by Sierra was to find a network composed of 
passive elements such that the output, f(t), was approximated by the 


normal curve 


Ur (a exp(-n°t“) -B<t<B QS), 


= 0) otherwise. 
Similarly, the error in this approximation at t = 8 1s 
© = exp(-n’B°). (1.4) 
The compression ratio, K, was defined as 
K = 20/28 . (1.5) 


Sierra assumed the same tolerable degree of error at both input and 


output pulse edges and equated (1.2) and (1.4) obtaining 
50 = iq (1.6) 


or E/n = IK 


The network transfer function was expressed as 
H(s) = F (s)/F.(s). (ia 
The two-sided Laplace transform of the input pulse was 


if exp(-e°t*) exp(-st) dt 


=_O 


F.(s) 


(yA/e) exp(s°/4e°). (1.8) 


For the output pulse, 


Fa(s) = (y#/n) exp(s@/4n°). (1.9) 

After substituting (1.8) and (1.9) into (1.7), the transfer function 
becomes 

2 
) 


he = iCenmennintlms) = (fe-)] soja) o (1.10) 


When (1.6) is substituted into (1.10), H{s) becomes 
_ gee 2 Ze 
H(s) = (1/K) exp[(i/K°) - 1] s/4e° . a, 


Since equation (1.11) still contains the parameter =, Sierra chose 


to normalize the input pulse by making € = 1, thus yielding 
H(s) = (1/K) exp[(1/K°) - 1] s*/4 . (1.12) 


At this point, Sierra introduced the parameter ¢, defining it as 


follows: 


(Ke -1)/4Ke Giaies 


a 
it 


It 


a ; Ko 12K (1.14) 


After substituting (1.13) into (1.12), H(s) reduces to 


H(s) = (1/K) exp(-9%s2) (1.15) 


As Sierra pointed out, the transfer function given by equation (1.15) 
is not realizable. Therefore, instead of the transfer function, another 
function was approximated. In this particular case, Sierra chose to 
approximate the magnitude function, noting the possibility of making 
this approximation and still obtaining accurate results in the time 
domain. 

Substituting s=jw into (1.15) and taking the magnitude, 


IH(jw)| = (1/K) exp(oew2) (i. Te) 
and ARG joe = O° . (1. 


The magnitude functions of the input and output pulses are 


JF.(jw)| = (Fr exp(-w°/4) (1.18) 
JF (jw)| =  (Wa/K) exp(-w°/4K°) (1.19) 
To obtain the frequency, w, , at which [est jl |F (Jw) | , Sierra 


equated (1.18) and (1.19), yielding the following: 
2 a poe Cee 
V7 exp(-w, /4) = (\m/K) exp(-w, /4K ) 


or w, = 2KEIn K/ (Ke - 17/2 (1.20) 


At this juncture, Sierra observed that the approximations to 
|H(jw)| must be bandwidth limited in order to obtain realizable transfer 
functions; that is, every approximation is valid only up to a frequency 
Was Sierra also pointed out that, since the network is to be realized 
with passive elements only, and with the same inpedance level for 
Wy > Ye » the necessity for introducing some attenuation factor, A, 
exists. 


Using equation (1.16), the magnitude-squared function can be 


obtained. It 1s 


HH(jw)|® = (1/K2) exp(2o%w*) (1.21) 


To approximate (1.21) by a ratio of polynomials in “, Sierra chose 
to use the continued-fraction expansion because of the very fast con- 


vergence it offers. 


Substituting z = 26%w* fai22)) 
into (1.21) yields 
il 2 
|H(z)|~ =. Gl Keexp(z) . (1.23) 


Expressing exp(z) as a continued fraction, the following is obtained: 














eZ) — 
THz 
2 
52 
Gia 
ez 
6 
| 
1 
nGEZ 2 
1 
Voz 
14 
ara (1.24) 


If N = number of terms taken in the expansion, the following approx- 


imations to exp(z) are obtained: 


POreNT= 4. t 
For N= 2: — 
° |-2 
7 2+Z 
BUD IE yas , etc, 


Thus, substituting for z yields the following approximations: 


For N= 1: |H(jw)|© = (1/K*) + 


1 - 26 w 
we 2 2, | + imi 
For N= 3 |H(G)| "= (ie a3 etc. 
1 - o> w 


Approximations to H(s)H(-s) were obtained by substituting s = ju. 


For N=1: H(s)H(-s) = (1/K) u 


For N= 2: H(s)H(-s) = (1/K°) 55 
1 + 26's 
Zee 
e 2, 1-05 
For N = 3: H(s)H(-s) = (1/K°) ——>-> 
ee 4°5° SeeLG. 


The poles and zeros for the approximations to H(s)H(-s) for N = 2 
through N = 9 are contained in Table I. The numerical values of these 
poles and zeros, that is, the values prior to being multiplied by 1/¢ 


will hereafter be referred to as the pseudo poles and Zeros. 
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ZEROS PUEES 


+-(0# 50. 70710678) 
141450) 1051) 
o o 
i-(#1.22474487+50) ;-(#0..3352200067:51 .054690705) 
1(+1.271229882 + §0.340625032 -(#0.340625032 +j1.271229882) 
{-(051.34867233) 


g #1 45534667 +j0.343560805 ) 
(+0.587391538 +j1.29829514) 


+(+1.52387182+50) (0+51.52387182) 


1(+1.47994076 +j0.5927200616) +(+0.5927200616 +j1.47994076) 
(+1.680548372+50) +(#0.79700116 +31.497202089) 
+(#1.64112335 +j0.59616025) +(20.250044024 +1.56720106) 
=(+1.6572801 +j0.801741003 ) +(+0.801741003 ol. 6SPAA01 ) 

+(+1.72038868 +j0.252045949) “-(#0.252045949 +j1.72038868) 


TABLE I. Poles and Zeros of H(s)H(-s) for N = 2 through 9. 


1] 


TT. COMPUTER SIMULATION 


In order to effectively analyze the Sierra Filter, several small- 
scale, special-purpose computer programs were required. Most of these 
programs were written to provide an output over a range of compression 
ratios for a given input. In addition to providing an optimum compres- 
sion ratio for a given transfer function, this procedure permitted the 
rapid evaluation of design parameters for a specified compression 
ratio. Computer programs were written to provide design information 
for all the Sierra transfer functions from N = 2 through N = 9 using 
both right-hand-plane (RHP) zeros and left-hand-plane (LHP) poles and 
LHP zeros and LHP poles. In some instances, therefore, minor modifica- 
tions to the basic program were required. These modifications, of 
course, depended upon the specific transfer function being considered. 
A brief description of the pertinent computer programs with the required 


modifications will now be delineated. 


A, PROGRAM NUMBER ONE 

This program generated sample values of the normalized input Gaussian 
pulse for various values of standard deviation. The mean was chosen to 
be zero for convenience. 

The output of this program provided the time, standard deviation, and 
value of the normalized pulse. A plot of the input pulse was also pro- 
vided. 

The primary purpose of this program was to provide the standard 
deviation for a specified pulse width and truncation error. This program, 
however, also provided sample values of the input pulse which could be 


used in spectrum analysis or convolution routines as required. 
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B. PROGRAM NUMBER TWO 

This program generated the poles and zeros for all the Sierra transfer 
functions listed in Table I for values of compression ratio varying from 
1.1 through 10.0 in steps of 0.1. In effect, this program performed the 
iterative multiplication suggested in columns 2 and 3 of Table I. The 


output, of course, provided the required poles and Zeros. 


C. PROGRAM NUMBER THREE 

This program computed and plotted the magnitude and phase of the 
Sierra transfer functions for values of compression ratio ranging from 
i.e through 3.0 in steps of 0.1. 

The inputs to this program were the pseudo poles and zeros listed in 
Table I and the desired compression ratios. The required iterations were 
performed by the program. 

The output of this program provided the compression ratio, the corre- 
sponding value of omega in radians and in degrees, and the corresponding 
magnitude and phase. A rapid means of comparing magnitude and phase for 
Specified compression ratios was made possible by this program output. 

Minor modifications to this program were required for each transfer 
function considered. These modifications related to the manner in which 
the transfer function was generated and will be included herein for con- 


venience. 


1. Modification for N = 2 


DO 70 I-1,100 
PST=2.*TK/SQRT(TK**2-1. ) 
$1=(0.0,0.70710678+W/PSI) 
S2=(0.0,-0.70710678+W/PSI1) 
PI=PSI*S1 

ieee SES 

W1=(1.0,0.0) 

W2=P1*P2 

T=(1./TK)*(W1/W2) 


Modification for N = 3 


DO 70 [=1,100 
PSIT=2.*TK/SQRT(TK**2-1. ) 
R1l=(-1.0,W/PST) 
$1=(0.0,1.0+W/PSI) 
S2=(0.0,-1.0+W/PSI) 

PSE SaeS 4 

ital obs Sy 

W1=Z] 

W2=P1*P2 
T=(1./TK)*(W1/W2 ) 


Modification for N = 4 


DO 70 f=1,100 

PS1=2,*TK/SORT (TK**2-1. ) 
R1=(-1.22474487 ,W/PSI) 
S1=(0.33522001,1.05469071+W/PSI) 
$2=(0.33522001,-1.05469071+W/PST) 
P1=PSI1*S] 

P2=PS1*S2 

P3=(1.0,0.0) 

P4=(1.0,0.0) 

Z1=PSI*RI 

Z2=(1.0,020) 

73=(1.0,0.0) 

74=(1.0,0.0) 

W1=Z1*Z2*73*Z4 

W2=P1*P2*p3*P4 

T=(1./TK)*(W1/W2) 


Modification for N = 5 


DO 70 I=1,100 

PST=2.. *1K/ SORT ( TK** 2-1 
R1=(-1.27122988 ,0.34062503+W/PSI ) 
R2=(-1.27122988 ,-0.34062503+W/PSI) 
S1=(0.34062503,1.27122988+W/PS1) 
S2=(0. 34062503 ,-1.27122988+W/PS1) 
Pl =p | 

od SD SY 

P3=(1.0,0.0) 

P4=(1.0,0.0) 

Z1=PSI*R} 

Z2=PSI*R2 

Z3=( 1050.0) 

Z4=(1.0,0.0) 

Wile 22" 25574 

W2=P1*P2*P3*P4 

T=(1./TK)*(W1/W2) 
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Modification for N = 6 


DO 70 I=1,100 

Bole 27 Ky SORMK** Zl) 

.. (-1.45534667 ,0.34356080+W/PST ) 
2=(-1.45534667,-0.34356080+W/PSI ) 
O20), 1); 34867233+W/PSI) 
2=(0.0,-1.34867233+W/PST) 
3=(0. 58739154, 1.2982951+W/PST) 

eae (0.58739154,-1.2982951+W/PSI) 

PSP oLsS | 

P2=P Sie SZ 

eee ea Sl FASS. 

P4=PSI*S4 

Z1=PSI*RI 

Z2=PSI*R2 

73=((0),.02 0) 

Z4=(1.0,0.0) 

WIl=Z1*Z2*73*Z4 

NCP Pers eS 

T=(1./TK)*(W1/W2 ) 


NO Chili aul © mea ea 


DO 70 [=1,100 

PST=2.*TK/SQRT(TK**2-1. ) 

a (-1.52387182,0.0+W/PSI ) 
2=(-1.47994076 ,0.59272006+W/PSI ) 
3=(-1.47994076 ,-0.59272006+W/PST) 
L=CORORIN 52387182+W/PSI) 
2=(0.0,-1.52387182+W/PST ) 
3=(0. 59272006, 1.47994076+W/PSI) 

rf (0. 59272006 ,-1.47994076+W/PSI ) 

Pr=PSies 

Beeb ols 52 

fo—PSI*S3 

P4=PSI*S4 

Z1=PSI*R1 

Z2=PSI*R2 

£o=Pon*R3 

Z4=(1.0,0.0) 

WI=Z1*2Z2*Z3*24 

W2=PI*P2*P3*P4 

T=(1./TK)*(W1/W2) 


Modification for N = 9 


DO 70 I[=1,100 

PSI=2.*1K/SORT(TK**2-1,, ) 

ie (-1.6572801,0.80174100+W/PSI) 
2=(-1.6572801,-0.80174100+W/PSI) 
3=(-1. 72038868, 0.25204595+W/PSI ) 
4=(-1.72038868, -0.25204595+W/PSI) 


lis 


S$1=(0.80174100,1.6572801+W/PSI ) 
S2=(0.80174100,-1.6572801+W/PSI) 
S3=(0.25204595,1.72038868+W/PSI ) 
S4=(0.25204595 ,-1.72038868+W/PSI ) 
Pl=PSI*S1 
P2=PSI*S2 
P3=PSI*S3 
P4=PS1*S4 
Z1=PSI*R1 
fZ=et RZ 
Z3=PSI*R3 
Z4=PSI*R4 
W1=Z1*Z2*73* 74 
W2=P1*P2*P3*P4 
T=(1./TK)*(W1/W2) 
D. PROGRAM NUMBER FOUR 

This program computed the maximum magnitudes of the Sierra transfer 
functions and the corresponding values of omega to five decimal places 
for values of compression ratio ranging from 1.2 through 3.0. 

The inputs to this program were the pseudo poles and zeros listed in 
Table I, the desired compression ratio and the values of delta, the step 
size for omega. This program computed the maximum magnitude as follows: 
Initially, the step size for omega was 0.1. When the maximum was passed, 
i1.e., omega (I) less than omega (I-1), the computer went to omega (I-2), 
changed the step size to 0.0] and continued the process. Any desired 
degree of accuracy could thus be obtained employing this technique. 

In many instances, it 1s necessary to multiply the original transfer 
function, H(s), by some constant in order to realize it. Specifically, 
in order to realize the constant-resistance lattice network, the require- 
ment that |H(jw)| be less than or equal to unity must be met. This 
program, therefore, provides an accurate value of the maximum magnitude 


so that these requirements can be fulfilled. 
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Since this program generated the transfer functions in the same 
manner as Program Number Three, the same basic modifications are appli- 


cable. These modifications will not be reiterated. 


E. PROGRAM NUMBER FIVE 

This program computes the numerator/denominator coefficients of the 
Sierra transfer functions for values of compression ratio ranging from 
1.] through 10.0 by performing the iterative multiplication of the fac- 
tored terms in the numerator and denominator. 

The input to this program was the pseudo poles and zeros. The output 
provided the required coefficients. This program was written to provide 
the coefficients using either LHP zeros or RHP zeros. Thus, the sign of 
the pseudo poles and zeros and the indicator value determined which 
coefficients were computed for a specified transfer function. 

Program Number Two generated the poles and zeros for the Sierra 


transfer functions and provided the transfer functions in the form 


(s+Z,)(S+Z,) ee ae 


= n 
Ms) = OVK) Tee US) SAT au 


Program Number Five, however, provided the transfer functions in the 


form 


a | 1s 6 
H(s) = (1/k) hat si 


Program Number Five was written, therefore, to provide additional 
design flexibility, since some network analysis programs can only accept 
transfer functions appearing in the form of (2.1) while others can only 


accept them in the form of (2.2). Also, analysis and synthesis frequently 


te? 


demand that both the poles and zeros and the coefficients of the transfer 
functions be available. 

Minor modifications to this program were required in order to calcu- 
late the numerator and denominator coefficients: for all the transfer 
functions. Since the same program was used to generate both the numera- 
tor and denominator coefficients, the FORMAT statements and the manner 
in which the coefficients were generated required modification. 

1. Modifications Required for Denominator Coefficients 

In order to label the output properly, the two FORMAT statements 
utilizing the Hollerith field must be changed when generating the deno- 
minator coefficients. 

2. Other Required Modifications 

By inspecting the poles and zeros contained in Table I, it is 
evident that the numerator and denominator coefficients for N = 9 and 
the denominator coefficients for N = 8 are generated in the same manner. 
These coefficients were generated as follows: 

70 =PSI=2.*1K/ SORMaI = 2— 1) 

Ms ((A**2 )+( Baggy) 

Z={ (C**2))1 (baercy) 
O=(PS1**4)*A1*A2 

lei RL Se (C*A1+A*A2) 
| 


PSI**2 )* (AMEE Aeon | 


53. U2 P Sie Aste) 


The numerator coefficients for N = 8 and N = 7 were generated in 
the same manner as shown in Program Number Five. 

The denominator coefficients for N = 7 and N = 6 were generated 
using the following modification: 

70 =PST=2.*TK/SQRT(TK**2-1. ) 
a icatamacam tales 28)) 
SO= (PSI**4 )*( (B**2)*Al ) 
SiRSs AAO) * 2a Oe BEA 2 ) 6 


Soar lee ) EB ta 
So=Bs 12,070 


The numerator coefficients for N = 6, the numerator and denomi- 
nator coefficients for N = 5, and the denominator coefficients for N = 4 
were generated as follows: 

70 PSI=2.*TK/SQRT(TK**2-1. ) 

Al=((A**2)+(B**2) ) 

SO=(PSI**2)*Al 

S1=PSI*2.0*A 

S$2=1.0 

$3=0.0 

The appropriate modification to this program, therefore, yielded 


the desired coefficients for the Sierra transfer functions. 


F. PROGRAM NUMBER SIX 

This program makes use of the DSL/360 (Digital Simulation Language) 
[Ref. 3] to provide the output for a given Gaussian input pulse and a 
given transfer function. The standard deviation of the input pulse can 
be varied several times in the same program, thus permitting a rapid 
comparison of the outputs versus input pulse width. Other parameters 
such as compression ratio and mean can also be varied as desired. 

The inputs to this program were the numerator and denominator 
coefficients of the transfer functions. These coefficients were generated 
by Program Number Five. A FORTRAN subprogram was written to generate the 
Gaussian input pulse with a truncation error of one percent. Program 
Number One, of course, provided the required standard deviation for a 


given pulse width. 


ihe 


III. SOME ASPECTS OF LATTICE SYNTHESIS 


In analyzing the Sierra Filter, synthesis was restricted to the 
symmetrical lattice. Similar synthesis procedures, however, could be 
applied to ladder development. In fact, in many cases, the symmetrical 
lattice can be decomposed into an equivalent ladder network [Ref. 4]. 
Figure 1 below shows the standard form of the symmetrical lattice and 


the bridge form, respectively. 





(a). Standard Form 


aes 


yO 


(b). Bridge Form 
Figure |. The Symmetrical Lattice 
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A. THE UNLOADED LATTICE 


Referring to Figure la, it is evident that 





mr fay t ope (3.1) 
Vay! 
Ts, Z,. = (2.3) 
1] I 4 
and Zs = Ye (3 4) 
Ke 2. theo 





Since the symmetrical lattice is being considered, Li 4200 and LioLoy: 


By applying elementary network analysis, it can be shown that 


cae 255 


NW 2"(ZaFPZb) 505) 


and i = i = 1/2 (Zb - Za) . (3.6) 


Consider now the transfer function, H(s), defined as follows: 


H(s) = Nl (3.7) 
> 7 Vays S I,=0 , 
or, from (3.1) and (3.2) 
z 
Hen = a (3.8) 
22 


Substituting (3.5) and (3.6) into (3.8) yields 


Hig) = BAe Fay, (2b: - 7a) (3.9) 
2 1/2 (Za + Zb) (ee Zay 


Dividing numerator and denominator of (3.9) by Zb, the transfer function 


can be written as 


2 | 


H(s) 


1 - Za/Zb 
1 + Za/Zb 


Solving for Za/Zb 


Za/Zb 
Thus, la = 
and Zb = 


1 - H(s) 
1 + H(s) 
1 - H(s) 
1 + H(s) 


(3) 10 


(Se 


2) 


Eel) 


are synthesis possibilities for the unloaded symmetrical lattice. 


B. THE CONSTANT-RESISTANCE LATTICE 


The symmetrical lattice terminated in a load of constant resistance 


R and possessing a driving point impedance of R at the input terminals 


will now be considered. 





Figure 2. 


Using Eqs. (3.1) and (3.2) and noting the fact that V 


with an input impedance of R. 


following equations are obtained: 


Vay) 


and -15R 


22 


Figure 2 below depicts this situation. 


Symmetrical lattice terminated in a resistance R 


= -I,R 


22! oR » the 


(3.14) 


(3. 1S 


Solving for L,, in (3.15) and substituting into (3.14) yields 





Vous = Lyly + Iyy (CAELY (3.16) 
Ti et S12 3 
fag 
Vay 
Solving (3.16) for i , the driving-point impedance is obtained. 
] 
V Jesu 
aj alee 
Zap S) 7 I, : Ly) Loy + — (3.17) 


Since Za, = Loo and Lyn = Loy > Eq. (3.17) can be written as follows: 
ne ee ule), 


Thus, for the driving-point impedance to equal R ohms, it 1S required 


that the right-hand portion of Eq. (3.18) equal R. In other words 


pp cea Ine RC 


17 (3.19) 


or Z = + R (320) 


Substituting (3.5) and (3.6) into (3.20) yielded 


1/4 (Za® + 27ZaZb + Zb°) = R° + 1/4 (Zb¢ - 2ZaZb + Za“) (3.21) 


or Tai = RC 


If R is normalized to 1 ohm, then Eq. (3.22) implies that, in order for 
the driving-point impedance to equal the constant value 1, Za and Zb 
must be reciprocal impedances. 

Using the normalized value for R and noting the fact that Vio! = -] 


the following relations were obtained from Eqs. (3.1) and (3.2): 
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V =e 


Va lala Zce 


and Von 2 2541, - LooVoo . (33 


Nea op) (oe 


aay 


24) 


Again, using the relations Ly = Loo and Lio = Lo) , the voltage 


transfer function defined by 


H(s) = ou 


was found. Multiplying (3.23) by Z and (3.24) by Z1, and 


12 
subtracting yielded 





2 a 
(29) = 2y9 + 24g) Voge = 2y0Vay) (3 
Ve Z 
“mmm (3 
L (24) = 22 + 2qy) 
. Substituting (3.5) and (3.6) into (3.27) gave 
Zb - Za 
ir = ory oS (3 
Since R is normalized to 1 ohm, Zb = 1/Za was substituted into Eq. 
obtaining 
_ | - Za 
WES)? eee 3 
pag Niereatliless) 
or Za = Toms y ie. 
and Zb = I/Za (3 


thus yielding a synthesis possibility for the constant-resistance 


symmetrical lattice. 
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a 45)) 


269) 


27) 


.28) 


(3525) 


oo 


30) 


ae), 


The restrictions that must be placed on H(s) in order for Za _ to be 
realizable will now be examined. 

First of all, H(s) can have no poles in the RHP. This restriction, 
of course, applies to all output-to-input transfer functions. Secondly, 
IH( Jw) | must be bounded and less than or equal to 1. A cogent discussion 
of the second requirement is given in Ref. 5. Thus, for Za and Zb 
to be positive real and, therefore, realizable, the following require- 
ments must be met: 

(1) H(jw) has no poles in the RHP or on the jw axis. 
ey may < 

Baycura [Ref. 6] outlines in detail a method for realizing the Sierra 
transfer functions as constant-resistance symmetrical lattice networks 
using the relations given in Eqs. (3.30) and (3.31). In fact, Ref. 6 
contains the realization of the Sierra transfer function for N = 9 and 
K = 2. Similar techniques will now be applied in realizing the Sierra 
transfer function for N = 8 and K = 2 using the relations given in Eqs. 


onlay and (P13). 
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IV. SYNTHESIS OF THE SIERRA TRANSFER FUNCTION FOR N=8 AND K=2 


Employing Eqs. (3.12) and (3.13) as a synthesis possibility for the 
unloaded symmetrical lattice, synthesis techniques for the Sierra trans- 
fer function for N = 8 and K = 2 were demonstrated. The outputs of the 
computer programs discussed in Chapter II were used extensively and 
greatly simplified the synthesis task. It is noteworthy that the syn- 
thesis techniques outlined herein can be readily applied to any of the 
Sierra transfer functions, the poles of which are located in the left- 


hand plane. 


A. SYNTHESIZING THE IMPEDANCE Za(s) 
The transfer function for N = 8 and K = 2 was divided by H(s) MAX 
to obtain 
3 2 
a,s ta,S ta,Sta 
ThE ds 1G MAX Bi eng (4.1) 
(s+ X1ls + X2)(S° Byes itay4)) 
The technique of resistance padding discussed by Guillemin [Ref. 7] 
could have been used, thus alleviating the necessity of dividing H(s) 
by H(s) MAX. This technique, however, can also be applied after the 
transfer function is normalized. 


Za(s), as given in Eq. (3.12) can now be written as 


Foals) 2=. = se ae (4.2) 
or Za(s) = Phas thse -_N(s) (4.3) 
a= ; 


Writing Za(s) as a ratio of polynomials 


26 


Ads* + M@s° + A2s¢ + Als + AO 


Za(s) 
Bs’ + B3s> + B2s¢ + Bls + BO 


(4.4) 
Program Number Two provided the poles and zeros for the original 

transfer function, H(s). Since the poles of Za(s) given in Eq. (4.4) 

are the same as the poles of the original transfer function, (4.4) was 


expanded in partial fractions as follows: 


2 2 
Ais. = SS 48 pe Seat OD : (4.5) 


5 foe t XZ 5 + X3s + X4 

Using the output of Program Number Two for N = 8 and K = 2, the 
factorization required by (4.5) was performed avd the values of X1, X2, 
X3, and X4 were thus obtained. 

Equations (4.4) and (4.5) are identities and must, therefore, be 
equal for all values of s. After equating (4.4) and (4.5) and clearing 
fractions, five equations in five unknowns were obtained. Program 
Number Five provided the denominator coefficients of Eq. (4.4), thereby 
Simplifying the task of clearing fractions. 

An existing subroutine was used to solve the above equations. The 


results were: 


A = O.47590151 = . 0.48 
h = I3.8173/@63 = £13.82 
C = -1.18232859 = -1.18 
» = -20.9721/061 = -20.9/7 
E = 1.70642708 = lee 


Substituting these values into Eq. (4.5) yielded 


2 2 
Za(s) = 0.48s" + 13.82 7 -1.18s~ - 20.97 1. 


rae Pe eee 
S ot S685) +015. 34 See}. 15Se+ $3.43 


2/ 


or Za(s) = Z,(s) 7 Z,(s) + Z..(s) 


(4.7) 


The equivalent impedance of a series LC circuit shunted by a resis- 


tance R was found to be 


Rs* + R/LC 


Z(s) 5 
Stee le Sot le 


(4.8) 


Thus, the first term in (4.6), Z,(s) , can be realized as a series LC 


circuit shunted by a resistance R. 


The second and third terms of (4.6) were written as follows: 


2 


2 
1,(s) + 2,(s) = -1.18s" ~ 20.97 + 1.71(s"#].15s#13.43) (4.9) 


ao Pe se oes 


Simplifying yielded 


Z 
Z(s) +Z,(s) = S795 
Sot lPos et alise4s 


Equation (4.10) was written as follows: 


2 
Z,(s) + Z,(s) - y oes + 2.03 + 
s + 1.158 + 13.43 
0.53s¢ + 2.03 
or Z,(s) 5 
S- + |, 5Setnlo las 
7 1.975 
and Z,(s) a. 7 


S$ taglloS at 13 ade 


tl aoe eca3 


(4.10) 


1.97s 
, (4.11) 


Ss + 1.18s +elamee 
(4.12) 


(4.13) 


Z,(s) is now in the form of Eq. (4.8) and can, therefore, be realized as 


a series LC circuit shunted by a resistance. 


The equivalent impedance of a parallel RLC circuit was found to be 
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Z(s) = —— (4.14) 
s- + (1/RC)s + (1/LC) 
Equation (4.13) has the same form as (4.14) and can, therefore, be 
realized as a parallel RLC circuit. Thus, Za(s) can be realized as 
three RLC circuits in series. Similar synthesis techniques can be 


applied to obtain a realization for Zb(s). 


Z5, 


V. CONCLUSIONS 


The Sierra Filters, like the linear phase filters discussed by Kuo 
[Ref. 8], produce increasingly complex transfer functions as the number 
of terms in the continued fraction expansion is increased. 

Several computer programs were written to provide a host of useful 
design parameters for virtually all the Sierra transfer functions ranging 
from N = 2 through N = 9 with the compression ratio and frequency varying 
over a Significant range. Although these programs do not provide all the 
design parameters that may be required, they frequently simplify the task 
of obtaining others. 

The use of small-scale programs that are designed to provide an out- 
put over a wide range of parameter variations can sometimes be more 
effective than the large-scale programs, some of which seem to handle 
parameter variations awkwardly and provide redundant information. 

Some synthesis techniques were delineated, primarily to demonstrate 
the value of these small-scale programs in synthesizing the transfer 
functions and, secondly, to demonstrate a synthesis technique that is 
relatively simple and can be applied to all the Sierra transfer functions 


that are realizable using passive elements. 
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10 
20 


30 
40 


90 


60 
70 


PROGRAM NUMBER ONE 
---TAYLOR---THESIS 


THIS PROGRAM GENERATES SAMPLE VALUES OF INPUT 
GAUSSIAN PULSE FOR VARIOUS VALUES OF STANDARD 
DEVIATION (SIGMA). THE MEAN IS EQUAL TO ZERO FOR 
CONVENTENCE. 


DIMENSION A(61),B(61),D(61,2) ,JXY(2) 

Z=] .00000000 

WRITE (6,20) 

FORMAT (1H1,4HTIME,7X ,5HSIGMA ,8X ,8HFUNCTION,// ) 
=-3.00000000 


CALCULATE THE INPUT PULSE VALUE AS A FUNCTION OF 
TIME AND ITERATE THE VALUE OF SIGMA. 

DO 40 I=1,61 

X=-(Z*T)**2 

F=EXP(X) 

A(1I)=F 

B(I)=T 

SIGMA=1./(SORT(2. )*Z) 


PRINT OUT TIME, STANDARD DEVIATION AND VALUE OF 
INPUT PULSE. 


WRITE (6,30) T,SIGMA,F 
FORMAT (1H ,F4.1,2F14.8) 
T=T+0. 100000000 

CONTINUE 


STORE VALUES IN A MANNER ACCEPTABLE TO THE PLOT 
ROUTINE. 


DO 50 I=1,61 
D(I,1)=B(1) 
D(1,2)=A(1) 
JXY(1)=1 
JXY(2)=2 


PLOT THE INPUT PULSE VERSUS TIME FOR VARIOUS VALUES 
OF STANDARD DEVIATION. 


CALL VPLOT (D,JXY,61,61,1,0,XMIN,XMAX, YMIN, YMAX ) 
Z=Z+0. 100000000 

IF(Z.EQ.4.0) GO TO 70 

GO TO 10 

STOP 

END 


3] 


FUNCTION F(X) 

F=EXP(X) "= 
RETURN 

END 
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PROGRAM NUMBER TWO 
---TAYLOR---THESIS 


THIS PROGRAM GENERATES THE POLES AND ZEROS FOR THE 
STERRA TRANSFER FUNCTIONS FOR VALUES OF COMPRESSION 
RATIO (TK) VARYING FROM 1.1 THROUGH 10.0 IN 
INCREMENTS OF 0.1. 


N=8 
COMPIBE XesRA PAZ in22., Sil 352 .RPIRZ 


WRITE (6,10) 
FORMAT (1H1,1X,2HTK,6X,3HPHI,14X,5HPOLES,19X, 


15HPOLES,19X,5HZEROS,19X ,5HZEROS,//) 


READ IN THE PSEUDO POLES AND ZEROS FOR THE 
TRANSFER FUNCTION. 


R1=(1.68054837 ,0.00000000) 
R2=(1.64112335,0.59616025 ) 
S1=(0.79700116,1.49720209) 
S2=(0.25004402,1.56720106) 
TK=1.1 


CALCULATE THE POLES AND ZEROS. 


DO 30 I=1,90 
PHI=2.*TK/SQRT(TK**2-1. ) 
P1=PHI*S1 

P2=PHI*S2 

Z1=PHI*RI 

Z2=PHI*R2 


PRINT QUT THE POLES AND ZEROS AND THE CORRESPONDING 
COMPRESSION RATIO. 


WRITE (6,20) TK,PHI,P1,P2,Z1,22 
FORMAT (1H ,F3.1,9F12.8) 
TK=TK+0. | 

CONTINUE 

STOP 

END 


Ss 
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PROGRAM NUMBER THREE 
~--TAYLOR---THESIS 


THIS PROGRAM COMPUTES AND PLOTS THE MAGNITUDE AND 
PHASE OF THE SIERRA TRANSFER FUNCTIONS FOR VALUES 
OF COMPRESSION RATIO VARYING FROM 1.2 THROUGH 3.0 
AND OMEGA VARYING FROM ZERO TO TEN RADIANS PER SEC. . 


COMPLEX T,R1,R2,R3 .R4 551552, c09n 1 Pc, h3,P4 Zl poem 
122 623,24 Wiewe 
DIMENSION A(100) ,B(100) ,C (100) ,D(100,2) ,E(100,2), 
1JX¥(2) 
10 READ(5,20) TK 
20 FORMAT (F10.9) 
N=8 
WRITE (6,30) N 
30 FORMAT (1H1,29HMAG AND PHASE VS OMEGA FOR N=,11,//1X, 
12HTK,6X , IOHOMEGA (RAD) ,5X, TOHOMEGA(DEG) ,7X, 
19HMAGNITUDE, 9X ,5HPHASE,// ) 
W=0.00000000 


GENERATE THE SIERRA TRANSFER FUNCTION IN A MANNER 
THAT PERMITS ITERATION OF OMEGA AND COMPRESSION RATIO. 


DO 70 I=1,100 

PS{=2.*1K/ SORMM TK=* elle) 
R1=CMPLX(-1.68054837 ,0.00000000+W/PST ) 
R2=CMPLX(-1.64112335 ,0.59616025+W/PST ) 
R3=CMPLX(-1.64112335,-0.59616025+W/PS1 ) 
S1=CMPLX(0.79700116,1.49720209+W/PSI ) 
S2=CMPLX (0. 79700116,-1.49720209+W/PSI) 
S3=CMPLX (0. 25004402 ,1.56720106+W/PS1I) 
S4=CMPLX (0. 25004402, -1.56720106+W/PS1) 
PHP Ses | 

P2=PoleSc 

Fo=PSiess 

P4=PSI*S4 

Z1=PSI*RI 

Z7=PSl RZ 

Z3=PSI*R3 

Z4=CMPLX(1.0;0.0) 

WIl=Z1*Z2*Z3*Z4 

W2=P1*P2*P3*P4 

T=(1./TK)*(W1/W2) 


COMPUTE THE MAGNITUDE OF THE TRANSFER FUNCTION. 


TMAG=CABS(T) 
A(1)=TMAG 
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100 


COMPUTE THE PHASE OF THE TRANSFER FUNCTION. 
PHASE=57 .29578* (ATAN2(AIMAG(W1 ) ,REAL (W1 ) ) 


3-ATAN2 (AIMAG(W2) ,REAL (W2) ) ) 


APPLY SCALE FACTOR TO PHASE. 
C(I )=PHASE/10.0 


OMEGA=57.29578*W 
B(I)=OMEGA 


PRINT QUT COMPRESSION RATIO, OMEGA IN RADIANS AND 
IN DEGREES, MAGNITUDE AND PHASE. 


WRITE(6,60) TK,W,OMEGA, TMAG, PHASE 
FORMAT(1H ,F3.1,4F16.8) 

W=W+0. | 

CONTINUE 

J=25 

ND IM=100 


PLOT THE MAGNITUDE AND PHASE VERSUS OMEGA. 


CONTINUE 

CALL VPLOT(E,JXY,0,NDIM,1,0,XMIN,XMAX, YMIN, YMAX ) 
[Fie EG. 3.0) GO TO 100 

GO TO 10 

STOP 

END 
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PROGRAM NUMBER FOUR 
~--TAYLOR---THESIS 


THIS PROGRAM COMPUTES THE MAXIMUM MAGNITUDE OF THE 
SITERRA TRANSFER FUNCTIONS FOR VALUES OF 
COMPRESSION RATIO RANGING FROM 1.2 THROUGH 3.0. 


COMPLEX T,R1,R2,R3,R4,S1,52,53,54,P1,P2,P3,P4,Z1, 
1Z2,23,Z24,W1,W2 
DIMENSION A(500) ,B(500),C(500) 
N=8 
WRITE (6,40) N 
40 FORMAT (1H1,27HHMAX VS OMEGA AND TK FOR N=,11,//2X; 
12HTK, 7X, ]OHOMEGA (RAD) ,6X,10HOMEGA(DEG) ,13X, 
1 8HMAG (MAX) ,/) 


READ IN DELTA, TK, AND J. DELTA IS THE STEP SIZE 
FOR OMEGA, TK IS THE COMPRESSION RATIO AND J IS AN 
INDICATOR. 


10 READ(5,20) DELTA,TK,J 

20 FORMAT(2F10.9,11) 

30 W=0.00000000 
B(1)=0.00000000 

90. «=«([=2 


GENERATE THE TRANSFER FUNCTIONS IN A MANNER THAT 
PERMITS ITERATION OF COMPRESSION RATIO AND OMEGA. 


60 PSI=(2.*TK)/SQRT((TK**2)-1.) 
R1i=CMPLX(-1.68054837 ,0.00000000+W/PST ) 
R2=CMPLX(-1.64112335,0.59616025+W/ PST) 
R3=CMPLX(-1.64112335,-0.59616025+W/PST ) 
S1=CMPLX(0.79700116,1.49720209+W/ PSI) 
S2=CMPLX(0.79700116,-1.49720209+W/ PSI) 
S3=CMPLX(0.25004402,1.56720106+W/PSI ) 
S4=CMPLX(0.25004402,-1 .56720106+W/PSI) 
PIEPsies| 
PZ=PS4- 52 
PoSpel.os 
P4=PSI*S4 
Z\=PSI*RI 
E7=P Slee 
Z3=PST*R3 
Z4=CMPLX(1.0,0.0) 

WI=Z1* 22% 274 
W2=P1*P2*P3*P4 
T=(1./TK)*(W1/W2 ) 


Io) 


MOMIOMmIOAND 


Mma 


70 


80 
90 


100 
110 


120 
130 


COMPUTE THE MAGNITUDE OF THE TRANSFER FUNCTION. 


TMAG=CABS(T) 
A(I)=W 
C(1)=57.29587*A(1) 
B(1)=TMAG 


COMPUTE THE MAXIMUM MAGNITUDE BY CONTINUALLY 
PASSING THROUGH IT AND REDUCING THE STEP SIZE 
FOR OMEGA. 


IF(B(I)-B(I-1)) 80,70,70 
W=W+DELTA 

[=I+] 

GO TO 60 

W=A(I-2) 

READ(5,90) DELTA, TK,J 
FORMAT (2F10.9,11) 

IF(J) 130,120,100 


PRINT OUT THE MAXIMUM MAGNITUDE, OMEGA IN DEGREES 
AND IN RADIANS AND THE COMPRESSION RATIO. 


WRITE(6,110) TK,A(I-1).C(I-1),B(I-1) 
FORMAT(1H ,F4.1,2F16.8,F22.8) 

GO TO 10 

[=I+] 

GO TQ 60 

STOP 

END 
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PROGRAM NUMBER FIVE 
---TAYLOR---THESIS 


THIS PROGRAM COMPUTES THE NUMERATOR/DENOMINATOR 
COEFFICIENTS OF THE SIERRA TRANSFER FUNCTIONS 
FOR VALUES OF COMPRESSION RATIO (TK) RANGING 
FROME Se MnCwGh lev0ean STEPS OF 0.1, 


READ IN THE PSEUDO POLES AND ZEROS, THE 
COEFFICIENT OF THE FOURTH POWER OF S AND THE 
INDICATOR, 1. 


10 READ (5,20) A,B,C,S4,1 
20 FORMAT(4F12.8,11) 
TK=1.1 
N=8 


IF I IS LESS THAN ZERO, RHP ZEROS ARE BEING 
USED Swe lie LSmEOUAL iC eZERO, EHP ZERGS 
ARE BEING USED. 


IF(I-1) 30,50,100 
30 WRITE(6,40) N 
40 FORMAT(1H1,28HNUMERATOR COEFFICIENTS USING,/1X, 
130HRHP ZEROS AND LHP POLES FOR N=,11,//2X,2HTK, 
19X,2HSO, 15X52HS1 5 5X .2HS2, one Sen OX ,ZHS 4a 


GENERATE THE COEFFICIENTS. 


70 PSI=(2.*TK)/SQRT(TK**2-1. ) 
All=NgrB* a2 y+ (0 saan) 
S0=(PSI**3)*A*AI 
S1=(PSI**2)*(A1+2.0*A*B) 
S2=PS1* (A+2.0*B) 

§3=1 .00000000 


PRINT OUT THE COEFFICIENTS AND CORRESPONDING 
VALUE OF COMPRESSION RATIO 
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WRITE(6,80) TK,SO0,S1,S2,S3,S4 
FORMAT(1H ,F4.1,5F16.8) 
IF(TK.GT.9.9) GO TO 10 
TK=TK+0. 1 

GO TO 70 

STOP 

END 
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STORAG 
TABLE 


PARAM 


CONTRL 
PRINT 
PREPAR 
GRAPH 
PRPLOT 
END 
PARAM 
END 
PARAM 
END 
PARAM 
END 
STOP 
FORTRAN 


PROGRAM NUMBER SIX 
---TAYLOR---THESIS 


THIS PROGRAM PROVIDES THE OUTPUT PULSE 
FOR A SPECIFIED INPUT AND TRANSFER FCN. 


IC (4) ,NUM(4) ,DEN(5) 
IC(1-4)=4*0,NUM(1-4)=0.5,-5.7305,... 
22.8391 ,-31.5524 ,DEN(1-5)=1,4.8360,... 
33. 02724167 iGS2, 200.0996 


TMU=3.0,SIGMA=0.9 
Y=GAUS (TMU, SIGMA, TIME) 
OUT=TRNFR(3,4,1C,NUM,DEN,Y) 
FINTIM=8,DELT=0. 1 
Orls UUT DEE 

Oe OUT 

TIME, OUT 

ONLY 


SIGMA=0. 70710678 
SIGMA=0 . 58925565 
SIGMA=0.41594517 


REAL FUNCTION GAUS(TMU,SIGMA,TIME) 
GAUS=EXP( (-0.5)*((TIME-TMU)/SIGMA)**2 ) 
IF(ABS(GAUS).LE.0.01) GO TO 5 

RETURN 

GAUS=0.0 

RETURN 

END 
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